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Abstract 

The dynamics of the quantum Fisher information of the parameters of the initial atomic state 
and atomic transition frequency is studied, in the framework of open quantum systems, for a 
static polarizable two-level atom coupled in the multipolar scheme to a bath of fluctuating vacuum 
electromagnetic fields without and with the presence of a reflecting boundary. Our results show 
that in the case without a boundary, the electromagnetic vacuum fluctuations always cause the 
quantum Fisher information of the initial parameters and thus the precision limit of parameter 
estimation to decrease. Remarkably, however, with the presence of a boundary, the quantum 
Fisher information becomes position and atomic polarization dependent, and as a result, it may be 
enhanced as compared to that in the case without a boundary and may even be shielded from the 
influence of the vacuum fluctuations in certain circumstances as if it were a closed system. 
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I. INTRODUCTION 


In estimation theory, the Cramer-Rao bonnd [1, 2] was proposed to describe how well one 
can estimate a parameter from the probability distribntion and the Fisher information is used 
to describe the precision limit. Since quantum mechanics is strongly related to probability 
theory in the sense that when we make quantum measurements on quantum mechanical 
systems, the observed outcomes follow a probability distribution, the Fisher information is 
readily extended to quantum regime and the inverse of the so-called quantum Fisher infor¬ 
mation (QFI) gives the lower bound of the error of the parameter estimation in quantum 
metrology [1-4]. With different models of the probe systems and different parameters to be 
estimated, many applications of quantum metrology have been done and some of those are 
of practical signihcance in quantum technology such as quantum frequency standards [5], op¬ 
timal quantum clock [6], measurement of gravity accelerations [7], clock synchronization [8], 
only to name a few. The central task in quantum metrology is to improve the precision of 
parameter estimation. Since a larger quantum Fisher information means better precision, 
increasing the QFI becomes a key issue in quantum metrology. A straightforward way, in 
addition to increasing the QFI, to enhance the precision when the probe systems are closed 
is by parallel measurements. Later it was shown that the use of correlated systems such as 
entangled states can also improve the precision of parameter estimation[5, 9-21]. On the 
other hand , however, interaction between a system and an environment is unavoidable in 
reality, and the quantum decoherence induced by such interactions may decrease the QFI 
and destroy the quantum entanglement in the probe system exploited to improve the pre¬ 
cision. In this regard. It has been shown that the interaction between a system and an 
environment usually makes the measurements noisy, which in turn degrades the estimation 
precision [22-41]. 

One environment which no system can be isolated from is the vacuum that fluctuates 
all the time in quantum sense. In the present paper, we are interested in how the vacuum 
fluctuations affect the QFI regarding the estimation of parameters of the initial states of 
a probe system which is modeled by a neutral polarizable two-level atom interacting with 
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electromagnetic vacuum fluctuations. We first examine how the decoherence caused by the 
vacuum fluctuations decreases the QFI as time evolves in an unbound space, then we ask what 
happens if the vacuum fluctuations are changed somehow, for example, by the presence of a 
reflecting boundary, i.e, can we protect the QFI from decreasing? We demonstrate that in an 
unbounded space, the QFI decreases exponentially with time. This behavior is similar to that 
observed in the amplitude-damping channel of decoherence derived in [40]. However, with the 
presence of a boundary, the situation changes dramatically in certain circumstances such that 
the QFI may even be protected from decreasing as if it were isolated from the environment. In 
other words, QFI may be shielded from quantum decoherence due to electromagnetic vacuum 
fluctuations in some special cases. At this point, let us note that the modihcation of vacuum 
fluctuations has been demonstrated to yield many interesting quantum phenomena such 
as the Casimir effect [42, 43], the light-cone fluctuations when gravity is quantized [44], the 
Brownian motion of test particles in an electromagnetic vacuum [45], the position-dependent 
spontaneous decay rates and geometric phase in an electromagnetic vacuum [46, 47]. 

II. QUANTUM FISHER INFORMATION AND DYNAMICAL EVOLUTION OF 
A TWO-LEVEL ATOM COUPLED WITH VACUUM FLUCTUATIONS 

For a given quantum state p{X) parametrized by an unknown parameter X, the unknown 
parameter can be inferred from a set of measurements, usually modelled mathematically 
by a set of positive-operator valued measures, on the state. After the optimization of the 
measurements and the estimator, a precision limit of the unknown parameter estimation is 
obtained [4] 

where N represents the repeated times and Fx denotes the quantum Fisher information of 
parameter X given by 

Fx = Tr(p(Jf)L2). (2) 

Here L is the so-called symmetric logarithmic derivative Hermitian operator, which satishes 
the equation dxp{X) = i{p(X),L} with {} standing for the anti-commutator. In the 
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orthonormal basis p{X) = t)e solved for and the quantum Fisher 

information can be written as [4] 


N 


Fx = 2Y1 


\{'lpm\dxP\^n)? 
Pm “1“ Pn 


( 3 ) 


For a two-level system, the state of the system can be expressed in the Bloch sphere repre¬ 
sentation as 

p= ^(J + n;-cr) , (4) 


where uj = (a;i, 072 , ^ 3 ) is the Bloch vector and cr = (ui, < 72 , < 73 ) denotes the Pauli matrices. 
As a result, Fx can be expressed in a simple form [40] 



u; < 1 , 


a; = 1 . 


( 5 ) 


Now let us hrst calculate the QFl for the initial state’s parameter estimation for an arbitrary 
state of the two-level atom 


IV’) = cos^|+)+e*'^sin^|-) , ( 6 ) 

where 9 and 0 correspond to the weight parameter and phase parameter, |-|-), |—) denote the 
excited state and ground state of the atom respectively, and the Bloch vector of the state 
can be represented as a; = (sin 6 ^ cos 0, sin 0 sin 0, cos 0). So the quantum Fisher information 
of 6 and 0 can be easily calculated as Fg = 1 and F^ = sin^ 9. Taking the atom as a closed 
system, whose evolution is governed by the Hamiltonian Hs = \ with uq denoting 

the transition frequency, one can easily show that the Bloch vector of the state with time r 
becomes u: = (sin 9 cos(0-l-a;or), sin 9 sin( 0 -|-a;o 7 ), cos 9), and the quantum Fisher information 
remains Fg = 1 and i0, = sin^ 9, which is time independent. Thus the unitary evolution does 
not change the quantum Fisher information of initial parameters. 

However, if the interaction with an environment is considered, the story may be different 
and the influence of the environment will in general be encoded in the atomic state with 
time. This is just what we are going to discuss next, and we consider how the QFl changes 
when the interaction with electromagnetic vacuum fluctuations is taken into account. For 
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this purpose, let us study a static polarizable two-level atom interacting with fluctuating 
electromagnetic flelds in vacuum, and in this case, the total Hamiltonian of the coupled 
system can be written as H = Hg + Hf + H', where Hg, which is given before, is the 
Hamiltonian of the atom. Hf denotes the Hamiltonian of the free electromagnetic held and 
its explicit expression is not required here. The Hamiltonian that describes the interaction 
between the atom and the electromagnetic held can be written in the multipolar coupling 
scheme as 

H\t) = -er ■ E(x(r)) , (7) 

where e is the electron electric charge, er is the atomic electric dipole moment, and E(a;) 
denotes the electric held strength. Let us note that, since both r(r) and E(a:) are not world 
vectors, the interaction Hamiltonian H' is ambiguous when we deal with atoms in motion. 
However, a manifestly coordinate invariant generalization of H' can be given [48]: 

H'(r) =-er^(r)F^^(x(r))M^(r) , (8) 

where is the held strength, r^(r) is a four-vector and its temporal component in the 
frame of the atom vanishes and its spatial components in the same frame are given by r(r), 
and is the four velocity of the atom. Since we have = (1,0, 0,0) in the frame of 

the atom, this extended interaction Hamiltonian reduces to that given by Eq. (7) in the 
reference frame of the atom. So we choose to work in this reference frame. Notice that we 
assume that the atom is static, as a result, the laboratory frame is equivalent to the frame 
of the atom. 

We let ptot = p(0) ® |0)(0| be the initial total density matrix of the system. Here p(0) 
is the initial reduced density matrix of the atom, which corresponds to the atomic state in 
Eqs. (6), and |0) is the vacuum state of the held. The evolution of the total density matrix 
Ptot in the proper time r reads as 

= (9) 

We assume that the interaction between the atom and held is weak. So, the evolution of the 
reduced density matrix p(r) can be written in the Kossakowski-Lindblad form [49, 50] 

^ = -^[i/=,,p(r)]+£b(r)], (10) 
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where 


1 V ^ 

^\p\ = 2^ Oij [2 ajp ai - CTjCrj p- p a^aj] . 
i,j=^ 

The coefficients of the Kossakowski matrix aij can be expressed as 

O'ij A6i^6j^ , 


( 11 ) 


( 12 ) 


with 

A = -[^(oio) + Gi—uJo)] , B = -[G{uJo) — Gi—uJo)] ■ (13) 

We define a two-point correlation function, G~^{x — x'), which is related to the two-point 
functions of the electromagnetic fields, {0\Ei{x)Ej{x')\0), as 

G+(i -x') = ^j2 {+l’'i|-)(-|rjl+) ( 0 |B,W£),(x')| 0 ) , (14) 

*J = 1 

and its Fourier and Hilbert transforms, ^(A) and /C(A), then follows 

/ oo p poo G(uj) 

dAr G+(At) , /C(A) = — / doo . (15) 

-OO J — oo ^ ^ 

By absorbing the Lamb shift term, the effective Hamiltonian Hqe can be written as 

Hes = ~ A5(a;o)]} <J^ , (16) 

where H is the effective level spacing of the atom. By applying Eq. (4) to Eq. (10), the Bloch 
vector with proper time r can be solved as: 

ui (r) = sin 6 cos(Hr -|- 0 ) , 

0)2 (r) = sin 6 * sin (Hr-I-0) , (17) 

'^ 3 ('r) = cos6e~^^'^ — ^(1 — . 


III. INFLUENCE OF VACUUM FLUCTUATIONS ON INITIAL PARAMETER 
ESTIMATION 

Let us now examine how the vacuum fluctuations affect the quantum Fisher information 
and thus the precision of the initial parameter estimation. Using the following electric two- 
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point function [51] 

he 

{0\Ei{x{T))Ej{x{T'))\0)o = - di&j) 

1 

(x — x'y + {y — y'Y + {z — z'Y' — {ct — ct' — isY ’ 
where £ ^ +0, 5' denotes the differentiation with respect to x' and the subscript 0 indicates 
the vacuum two-point functions in the unbounded space. Applying the trajectory of the 
atom 


((t) = T , x(t) = Xo , y{T) = ya , z(t) = Zo , 


(19) 


we find that the electric-held two-point functions can be written as 


{0\Ei{x{T))Ej{x{T'))\0)o = 


he 


6. 




tt^Sq (c At — iey ’ 

So the Fourier transform of the correlation functions can be calculated as 


6'“’(A) = 5j 


Sireohc^ 


9(A), 


( 20 ) 


( 21 ) 


with 9{X) being the standard step function. Let us note that here we let £ = 0 after the 
calculation. The coefficients of the Kossakowski matrix Ojj and the effective level spacing of 
the atom are now given by 

yl(o) ^ ^(0) ^ 


= 0)0 + 


7o 

2ttuJq 


P 




. ho 

4 ’ 

(22) 

f 1 M 

(23) 

\uj + ojq 00 -OOoJ 


Here 70 = e^|(—|r|-|-)pa;o/37r£o^c^ denotes the spontaneous emission rate in vacuum without 
boundaries, As a result, the quantum Fisher information of the initial weight and phase 
parameter become 

(24) 


E^P'^ = sin^ 9 e 


- 70 T 


and 

^ g-7or _ ^25) 

This shows that the QFI of both weight and phase parameters decreases exponentially with 
time due to the decoherence caused by the interaction between the atom and the fluctuating 
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vacuum and the decay rate equals the spontaneous emission rate of the atom in vacuum. 
Therefore, when the measurement time is larger than the relaxation time I/ 70 , the precision 
of the estimation is greatly damaged. This kind behavior of the QFI is the same as that in 
the amplitude-damping channel described in [40] as expected. 

Since vacuum fluctuation will be modihed if we set a boundary in the vacuum, we may 
wonder how the presence of a boundary which confines the vacuum fluctuations of the held 
affects the quantum Fisher information of the initial parameters with time. Now we consider 
the case with the presence of a rehecting boundary ai z = 0 The electric two-point 
functions in this case can be expressed as a sum of the Minkowski vacuum term and a 
correction term due to the boundary: 



where the {Ei{x{T))Ej{x{r')))o is the two-point function in the unbounded vacuum which 
has already been calculated above and 



gives the correction due to the boundary. Here n = (0, 0,1) is the unit vector normal to the 
boundary. By applying the trajectory in Eq. (19), the boundary term of the electric-held 
two-point functions can be written as 


(0|E^(x(r))E^(x(r'))|0)b = (0|Ey(a:(r))Ey(x(r'))|0)b 


he c^At^ + 4^2 


[ (c Ar — iey — 4 :Zq]^ ’ 


(0|E^(a;(r))E^(a;(r'))|0)b 


he 


1 


(28) 


^ ^ TT^eo [ (c Ar — ie)2 — 4^0]^ 

The Fourier transform of the correlation functions can be calculated as [47]: 



0 


(29) 


^ Let us note here that external forces may be needed to balance the Casimir-Polder force to keep a polar¬ 
izable atom at a fixed position near a reflecting boundary. 








with 


3c3 


fx{X,zo) = fy{X,zo) = 


2,Xzn 


- COS ■ 


2Xz, 


3c3 


UX,Zo) - g^3 3 


2Xzr] 


■ COS 


c 

2Xz() 


+ 


AX'^Zq \ . 2Xzo 
—-zz — 1 sin- 


2Xzr\ 


Sin 


(30) 


(31) 


Thus the coefficients of the Kossakowski matrix a^- and the effective level spacing of the 
atom now become 

A = B = ^{l-Y,»^M^o,Zo)), (32) 


where cij = |(—|rj|+)p/|(—|r|+)p. Physically, a* represents the relative polarizability and 
they satisfy = 1- We let f{uo,Zo) = '^^aifi{uo, Zq). As a result, the quantum Fisher 

information of the initial weight and phase parameters can be expressed as follows 


= sin^ e 


(34) 


and 


Pg = g-7o(l-/(KiO,2o))r 


(36) 


So, the decay rate of the QFI is now modified by the factor, 1 — f{u)o,zo), as compared to 
the unbounded case. Let us now first examine the asymptotic behaviors of the QFI, that 
is when the atom is placed very far from the boundary or very close to it. When zq ^ oo, 
fi{uo,zo) 0, the unbounded case is recovered as expected. When zq —?■ 0, fx{zJo,zo) = 
fy{uJo,Zo) = —fz{zoo,Zo) = 1, so atoms with different polarization will behave differently. 
When the polarization is along the 2;-axis, i.e., ax = cty = 0, the decay rate becomes double 
of that in the unbounded Minkowski vacuum case, which makes the QFI decay even faster 
than that without the boundary. However, when the polarization is in the xy plane, i.e., 
CK^ = 0, the decay rate becomes zero, which means that the QFI is totally protected from 
electromagnetic vacuum fluctuations for transversely polarizable atoms extremely close to 
the boundary as if it were isolated. For an isotropic polarization, cKj, = ck^ = ck^ = |, the 
decay rate is |7o, so the QFI decreases slower and the precision is enhanced by the presence 
of the boundary as compared to the unbounded case. 
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For a generic position and polarization, the QFI may be decreased, or enhanced as com¬ 
pared to the unbounded case. This can be seen from the fact that fx decays in an oscillatory 
manner from 1 to 0 with increasing atomic distance, while fz does from —1 to 0. In general, 
there exist some special positions where ^o) = 0. The quantum Fisher informa¬ 

tion of the initial parameters of the atom at these positions takes the same form as that 
in the unbounded case and the boundary effects disappear. In regions between these spe¬ 
cial positions, we have either ^o) > 0 or < 0. As a result, the quantum Fisher 

information is either enhanced or decreased as compared to that in the unbounded space. 
To show the properties we described above graphically, we plot, in Fig. (1), 1 — f{u)o,Zo) 
as a function of zo for cx = (1, 0, 0), (0, 0,1), (1/3,1/3,1/3), which corresponds to the paral¬ 
lel, vertical, isotropic polarization cases respectively. The oscillatory behavior displayed in 




Zo 


Zo 


Zo 


FIG. 1: 1 — / as a function of zq for a = (1,0,0), (0,0,1), (1/3,1/3,1/3) respectively. Here zq is 

in the unit of c/ujq. 

Fig. (1) is related to the stationary modes of the fields as a result of the superposition of 
the propagating incident and reflected modes. The oscillations of both the tangential and 
normal components of the field lead to the different oscillatory behaviors of the decay rates 
of both the atoms polarizable in the directions parallel and vertical to the planar boundary. 
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IV. EFFECTS OF VACUUM FLUCTUATIONS ON ATOMIC FREQUENCY ES¬ 
TIMATION 


As we have demonstrated, the presence of a reflecting boundary may protect the esti¬ 
mation precision of the initial parameters from the influence of the environment in certain 
circumstances. We may also wonder how it affects other parameters’ estimation such as the 
frequency estimation. Now, our unknown parameter becomes the atomic frequency Uq. By 
using Eq. (5) and Eq. (17), the quantum Fisher information of Uq in the unbounded case 
and the bounded case can be calculated, after dropping the terms which are of higher order 
in terms of the hne structure constant in the effective level spacing of the atom, to the hrst 
order as 


F® ^sin^0e"'^°W^ + 16, ^ 

° ' ou, 


d-fo 


1 + 


sin^^ 


+ cos 9 (cos 9 + 2) e 


g-70-^2 


(36) 


and 


^ sin^ 9 e-T'oAi-/(‘^o,^o))^2 ^ 


16 


^[7o(l-/(^o,^o))] ]' 
duo f 


X 


1 + 


sin^ 9 


+ COS0 (cos0 + 2) e-^oAi-/(-oxo)) 


X 


g-7or(l-/(wo,^o)).j-2 


(37) 


So, the maximum quantum Fisher information is obtained when 9 = which is given, after 
dropping the higher order term, by 

FiO) ^ g-7or^2 


and 




(39) 


Applying the equation = 0, we hnd that the maxima of Fi^^ and ^ and 

— 77-726“^, can be reached at = — and r = —r. —J -vr- When zq —)■ 00, 

7o[1-/AoXo)]^ ’ 70 7o(l-/(too,2o)) ^ ’ 

fi{uo, Zq) —)■ 0, the bounded vacuum case reduces to the unbounded case as expected. When 
Zq —!■ 0, /x(i^o,^o) = fyi^OyZo) = —fz{(^o,zo) = 1, thus different polarization directions lead 
to different results. When the polarization direction is along the 2; axis, i.e., ax = cty = 0, 
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the optimal measurement time and the maximal quantum Fisher information in the bounded 
vacuum case become 't = ^ and The optimal measurement time is half of that 

in the unbounded vacuum case and the maximal quantum Fisher information is a quarter 
of that in the unbounded case. When the polarization direction is in the xy plane, i.e., 
= 0, both the optimal measurement time and the maximal quantum Fisher informa¬ 
tion in the bounded case approach to infinity, and the quantum Fisher information with 
time r approaches to r^, which is consistent with the quantum Fisher information in uni¬ 
tary evolution ijJ = (cos(0 -|- cuot), sin(0 -|- uqt), 0). When polarization direction is isotropic, 
ctx = cty = <^z = the optimal measurement time and the maximal quantum Fisher in¬ 
formation in the bounded vacuum case become r = — and Fb„ = which are both 

larger than those in unbounded vacuum case. For an arbitrary polarization, in regions where 
/(cuq, ^o) > 0, both the optimal measurement time and the maximal quantum Fisher infor¬ 
mation are enhanced as compared to the unbounded vacuum case, while in regions where 
f{uo,zo) < 0, they are depressed. 

Let ns note that in estimation of a transition frequency, the total time T of probing is also 
a resource. For a given T S> we can improve the precision by accomplishing a sequence 
of measurements [23]. We let iV = T denote the repeated times of the measurements. Then, 
according to Eq. (1), the uncertainty of the atomic frequency satisfies 


lAwol > 


1 

/VF5 


1 


1 

-y/2"7-e-707-(l-/(aJo,2:o)) 


(40) 


Applying = 0, we obtain the minimum uncertainty of Uq, i.e., \Auo\min = 

^ e'ro(i-f^uio,zo)) ^ optimal sequence measurement time r = . As a result, 

the optimal repeated times of measurements in a given total time of probing to maximize 
the precision is T7o(l — /(cuq, ^o))- In regions where /(cuq, Zq) > 0, the precision is enhanced 
and the measurement times we need to obtain the maximum precision is less than that in 
the unbounded case. In regions where /{ojq, zq) < 0, the precision is degraded and the mea¬ 
surement times become more than that in the unbounded case. It could also be possible to 
further improve this precision by the use of an entangled probe system following [23], the 
detailed analysis of which is left as future work. Here, we give a very brief comment on the 
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issue. Take N maximally entangled atoms in vacuum for example and let ptot = p(0) <8 |0)(0| 
be the initial total density matrix of the system, where p(0) is the initial reduced density 
matrix of the N maximally entangled atoms. Then we can use Eq. (9) to study the evolution 
of the total state. After obtaining the state of the atoms in time, we can in principle calculate 
the quantum Fisher information using Eq. (3). Because of the indirect interaction between 
atoms caused by the held, it is not easy to hnd an analytical result of the reduced density 
matrix of the N atoms in time r. But in the special case when all the coefficients of the 
dissipative part in the evolution equation of the reduced density matrix vanish, all the atoms 
behave like closed systems and the Heisenberg precision limit is protected from deterioration. 
This condition can be approximately fulhlled if atoms are transversely polarizable and very 
close to the boundary. 


V. CONCLUSION 

In summary, we studied the dynamics of the quantum Fisher information for the atomic 
parameter estimation for a static polarizable two-level atom coupled in the multipolar scheme 
to a bath of huctuating vacuum electromagnetic helds without and with the presence of a 
rehecting boundary. When we estimate the parameters of initial atomic state, we hnd that, 
in the case without a boundary, the electromagnetic vacuum huctuations always cause the 
quantum Fisher information of the initial parameters and thus the precision limit of param¬ 
eter estimation to decrease. However, with the presence of a boundary, the quantum Fisher 
information becomes position and atomic polarization dependent, and thus the precision of 
the initial parameter estimation may be decreased, enhanced or remain unchanged as com¬ 
pared to the case without a boundary depending on the position and polarization . When 
the atom is extremely close to the boundary and is transversely polarizable, the quantum 
Fisher information may even be shielded from the influence of the vacuum fluctuations and 
remains constant with time as if it were a closed system. For the estimation of the atomic 
frequency, there exist a maximum quantum Fisher information and optimal measurement 
time, which can also both be enhanced or decreased as compared to the case without a 
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boundary. 
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